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There has been much interest in so-called SIC-POVMs: rank 1 symmetric in- 
formationally complete positive operator valued measures. In this paper we 
discuss the larger class of POVMs which are symmetric and informationally 
complete but not necessarily rank 1. This class of POVMs is of some indepen- 
dent interest. In particular it includes a POVM which is closely related to the 
discrete Wigner function. However, it is interesting mainly because of the light 
it casts on the problem of constructing rank 1 symmetric informationally com- 
plete POVMs. In this connection we derive an extremal condition alternative 
to the one derived by Renes et al. 
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1. Introduction 

There has been much interest in rank 1 symmetric, informationally complete 
positive operator valued measures [IlEllSlllligEllTllHlElIIOliniinillS]: SIC- 
POVMs, as they are often called. In d-dimensional Hilbert space these are systems 
of <P operators Er = {^/d)Pr such that each Pr is a rank 1 projector and 

TT{PrPs)^J^^{l+6rs) (1) 

for all r, s. In that case it can be shown that X]r=i -^'^ ~ operators Er 

constitute a POVM. Moreover the POVM is informationally complete P^fT^fTCl 
[T71 UHl nil HO] (meaning that an arbitrary density matrix p is completely specified 
by the probabilities Tv{Erp))- POVMs of this kind have been constructed [Ullllll 
m EH [22] (analytically and/or numerically) for every dimension d < 45. It is still 
an open question whether they exist in dimensions > 45. 

The purpose of this paper is to discuss POVMs which are still symmetric, in the 
sense that 

i:T{ErEs)=a + I35rs (2) 

for fixed numbers a, and informationally complete, but which are not assumed 
to be rank 1. We will refer to such POVMs as SI-POVMs ("S" for "symmetric", 
"I" for "informationally complete"). SI-POVMs which are also rank 1 we will refer 
to as SI(l)-POVMs (so an SI(1)-P0VM is what in the literature is often called a 
SIC-POVM). 

SI-POVMs are of some independent interest. In particular, we will show in 
Section [7] that the discrete Wigner function is closely related to a POVM of this 
type. However, our main reason for studying them is to gain additional insight 
into the problem of constructing SI(l)-POVMs. To that end we derive an extremal 
condition alternative to the one used by Renes et al in their numerical work. 

The plan of the paper is as follows. In Section [2] we discuss some geometrical 
features of quantum state space which will be needed in the sequel. In Section [3] 
we relate this discussion to the problem of devising a tomographical procedure 
which is, in some suitably defined sense, optimal. In Section U] we prove a theorem 
characterising the structure of an arbitrary SI-POVM. In Section [5] we specialise to 
the case of SI-POVMs covariant under the Weyl-Heisenberg group (or generalized 
Pauli group as it is often called). We show that such POVMs have a very simple 
representation in terms of the Weyl-Heisenberg displacement operators. In SectionlH] 
we turn to the problem of constructing SI(l)-POVMs, and derive an extremal 
condition alternative to the one derived by Renes et al [T]. Finally, in Section [7] we 
construct an SI-POVM which is closely related to the discrete Wigner function. 



2. The Block Body 

Let Ti he a, d dimensional Hilbert space, and let V be the space of density 
matrices defined on 7i. If d = 2 it is well known that V can be identified with the 
Bloch sphere. To be specific: let B be the imit ball in R-^ (i.e. the set of vectors 
G 'M? having length < 1). Then a 2 x 2 complex matrix p is a density matrix if and 
only if it can be written in the form 

p=\{l + h.cr) (3) 

where h € B (the Bloch ball) and cfi, <j2, ^3 are the Pauli matrices. 

With the appropriate modifications this construction can be generalized to higher 
dimensions [HIlSlIMlIlSIlSIllIiaiiaiSni- Let su(d) be the - 1 dimensional 
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real vector space consisting of all trace zero Hermitian d x d complex matrice^. 
Let B be the convex subset consiting of all B e su(d) for which B > —1. Then a 
d X d matrix p is a density matrix if and only if 

/5-^(l + B) (4) 

for some B G B. We refer to B as the Bloch body, and to its elements as Bloch 
vector^ 

It is convenient to define an inner product on su(d) by 

iB^,B,) = ^^TriB.B2) (5) 

for all -81,-82 G su(o?) (so {Bi, B2) is just the Hilbert-Schmidt inner product rescaled 
by the factor -jjj^)- Let 



l|B|| = V(^'^) (6) 

be the corresponding norm. 

If d = 2 a vector B G sn{d) is a Bloch vector if and only if ||-B|| < 1. Moreover, 
the corresponding density matrix is a pure state if and only if ||-B|| — 1. For d > 2 
the situation is more complicated. Let Bi and Bo be the balls 



Bi 



{Besu{d):\\B\\<^} (7) 

50 = {Besu(d): IIBII < 1} (8) 

and let 

51 = {sesu(d): 11511 = ^} (9) 

5o = {s e su((i): ||B|| = 1} (10) 
be the bounding spheres. Then 

BiQBCBo (11) 

It can further be shown [531 HSl 123 that Bi and Bo are respectively the largest and 
smallest balls centred on the origin for which this is true. Specifically: 

(1) If r > l/(d - 1) there exists B e su{d) such that ||i?|| = r and B B. 

(2) If < r < 1 there exists B G B such that ||-B|| = r. 

Moreover a Bloch vector B G B corresponds to a pure state if and only if it has 
norm = l{i.e. if and only if it G i3 n So)- 

It is worth noting that Bengtsson and Ericsson [B] have proved a stronger result: 
in any dimension for which either a full set of MUBs (mutually unbiased bases) or 
an SI(1)-P0VM exist Bi is the largest ellipsoid which can be inscribed in B. 

If d — 2 we have Bi = B = Bo and B Ci So — So, so the Bloch body has a very 
simple geometrical structure (it is just a ball of radius 1 centred on the origin, with 
the pure states comprising the boundary). For d > 2 these relations no longer 
hold, and the geometry is much harder to appreciate intuitively. One gets some 
additional intuitive feeling for the geometry, at least in low dimension, by looking 



^su(ci) is the Lie algebra for the special unitary group SU(d). This group theoretical fact is 
highly relevant to the problem of characterizing the geometry of quantum state space 24 , 25] I2til 
1271 1281 [29I I30j . However, it will play no part in the considerations of this paper. 

^What we are calling Bloch vectors are of course matrices. Some authors introduce a standard 
basis for s\i{d) at this point and reserve the term "Bloch vector" for the components of B in 
that basis (as has been the long-standing practice in the 2 dimensional case — see Eq. ^ above). 
However, it appears to us that this makes the notation needlessly complicated. 
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at the 2-dimensional sections of B which have been calculated [531 113 Ull HI] for 
d = 3 and 4. 

Let B be any vector e So (not necessarily a Bloch vector). An immediate 
consequence of Eq. pT|) is that xB e B whenever \x\ < — 1). Kimura and 
Kossakowski pSj have proved some much stronger results. In the first place they 
have shown 

Theorem 1. Let B be any vector £ So (not necessarily a Bloch vector). Let — to_ 
be the smallest eigenvalue of B and let m-|_ be the largest (so — r7i_ < B < ). 
Then 

(1) The quantities m± satisfy the inequalities 

1 <m- <d-l (12) 

and 

1 <m+ <d-l (13) 

Moreover to_ = I if and only if m+ = d — 1, and TO-|_ = 1 if and only if 
ru- = d ^ 1 

(2) B is a Bloch vector (in fact the Bloch vector corresponding to a pure state) 
if and only if m- = 1. Similarly —B is a Bloch vector (in fact the Bloch 
vector corresponding to a pure state ) if and only if m+ = 1 . 

Proof. See Kimura and Kossakowski [29] . □ 

Theorem [T] characterizes the vectors G BOSo {i.e. the Bloch vectors correspond- 
ing to pure states) in terms of their eigenvalues. The next theorem relates the 
diameter of the Bloch body in the direction B to the eigenvalues of 13. 

Theorem 2. Let B and m± be as in the statement of Theorem]^ and let x G R. 
Then xB & B if and only if 

< a; < (14) 

Proof. See Kimura and Kossakowski [35] . □ 

Remark. As Kimura and Kossakowski point out, it follows from Theorems [1] and [2] 
that a point where the boundary of B touches the outer sphere So is always dia- 
metrically opposite to a point where the boundary of B touches the inner sphere Si 
(and conversely). 

We conclude this section by proving a theorem which shows that, instead of 
considering the eigenvalues (as in Theorem [l]), one can use the quantity Tr(B3) to 
tell whether a vector B S 5o is the Bloch vector corresponding to a pure state. We 
first need to prove 

Lemma 3. Let P be any dxd Hermitian matrix (not necessarily a positive matrix). 
Suppose 

Tr(p2) = 1 (15) 

Then 

Tr(p3) < 1 (16) 
with equality if and only if P is a one dimensional projector. 
Remark. It is not assumed that Tr(P) = 1. 
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Proof. Let Ai, A2, . . . , Ad be the eigenvalues of P (not necessarily distinct). In view 
of Eq. (HSl) 

Define 

d 

« = E|A.|3 (18) 

It follows from Eq. ^T7\ that \Xr\ < 1 for all r, and consequently that 1 — |Ar| > 
for all r. So 



1-« = E(A^-|A.|3) 



r=l 
d 

= Ea^(1-|A.|) 

r=l 

> (19) 

with equality if and only if A^(l — |A,. |) = for all r. Consequently 

K < 1 (20) 

with equality if and only if 

A^(1-|A.|)=0 (21) 

for all r. 

It is now immediate that 

Tr(p3) < |Tr(p3)| < k < 1 (22) 

Suppose 

Tr(p3) = 1 (23) 
Then it follows from Eq. (1^^ that k — 1 which means, in view of Eqs. ^U\i and (PTjl . 
that A^(l — |A,.|) = for all r. Consequently, for each r, \Xr\ =0 or 1. The fact that 
J2r -^r ~ 1 then implies that |Ar| — 1 for exactly one value of rand — for all the 
others. Since, by assumption, A^ = 1 we must actually have A^ = 1 for exactly 
one value of r and = for all the others — implying that P is a one dimensional 
projector. 

If, on the other hand, P is a one dimensional projector it is immediate that 
Tr(p3) = 1. □ 

We are now in a position to prove our main result: 

Theorem 4. Let B be any vector G So (not necessarily a Block vector). Then 

(1) The quantity Tr(B^) satisfies the inequalities 

- d{d - l)(d - 2) < Tr(p3) < d{d - l)(d - 2) (24) 

(2) The upper bound in Inequalities {24-\j is achieved if and only if B (z B H So 
(and is therefore the Bloch vector corresponding to a pure state). 

(3) The lower bound in Inequalities \24-\j is achieved if and only if ^B ^ BOSo 
(and is therefore the Bloch vector corresponding to a pure state). 

Proof. The fact that 13 £ So means 

Tr(p2) = d{d - 1) (25) 

Define 

P± = i(l±P) (26) 
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Then Eq. (gS]) implies 

Tr(P|)-^(d + Tr(B2)) =1 (27) 

We may therefore use Lemma [3] to deduce 

^(fi + 3Tr(B2)±Tr(S^)) = Tr(P|) < 1 (28) 

with equahty if and only if P± is a one dimensional projector. In view of Eq. (|25p 
this means 

Tr(S3) < 2) (29) 

with equality if and only if P+ is a one dimensional projector, and 

Tt{B^) > -d{d - l){d - 2) (30) 

with equality if and only if P_ is a one dimensional projector. But is a one 
dimensional projector if and only if _B G B H Sq, and P_ is a one dimensional 
projector if and only if —B (z B D Sq. The claim is now immediate. □ 

3. Bloch Geometry and Tomography 

The geometry of the Bloch body is intimately related to the problem of devising 
measurement schemes which are, in some suitably defined sense, tomographically 
optimal. The connection works both ways. On the one hand knowledge of the 
geometry tells us what measurement schemes are possible. On the other hand a 
knowledge of possible measurement schemes provides important insight into the 
geometry. In this section we summarize the Bloch geometrical aspects of two such 
measurement schemes: namely, schemes based on a full set of mutually unbiased 
bases or MUBs [3ll[3li[ll[3Tl[13[S3[31[13[li[n3[li[iaiin] and schemes 
based on SI(l)-POVMs [HiaEllllElElITlIHlllllinillllinillS! (or SIC-POVMs as 
they are often called). Much of the material in this section amounts to a review 
of the relevant parts of Bengtsson and Bengtsson and Ericsson |6j , but using a 
slightly different terminology and notation. 

We begin with the case of a full set of MUBs. Suppose one has a large number of 
copies of a d-dimensional quantum system, all presumed to be in the same quantum 
state. Suppose one takes a fixed von Neumann measurement having d distinct out- 
comes, and performs it on many copies of the system. Suppose one then identifies 
the relative frequencies obtained with the corresponding probabilities. This will 
give one d probabilities pi,p2, ■ ■ ■ ,Pd- Taking into account the normalisation con- 
dition X]r=i^'^ ~ ^ ^^^^ means one has d — 1 independent numbers. On the other 
hand a full specification of the quantum state requires d^ — I independent numbers. 
It follows that if one wants to perform tomography using only von Neumann mea- 
surements one needs to divide the set of copies of the system into a minimum of 
d + 1 subsets, and to perform different von Neumann measurements on the copies 
belonging to different subsets. We will refer to a measurement scheme based on 
the minimum number of d -I- 1 different von Neumann measurements, each having 
d distinct outcomes, as a minimal Von Neumann scheme. 

The question now arises: what is the best way of choosing the d + 1 different 
measurements in a minimal von Neumann scheme? Let Pf,PJ, . . . be the d 
orthogonal, one dimensional projectors describing the r^^ measurement and let 
Bl,B2, ■ ■ ■ be the corresponding Bloch vectors. So 

p:^1{i + b:) (31) 

for all a, r. Notice that, whereas in Section [5] we used Bloch vectors to describe 
quantum states, now we are using them to describe quantum measurements. Notice 



6 



also that the fact that the are all one dimensional projectors means that the 
vectors IB^ all He on BCiSo- 
The orthonormality condition 

Tt{P:PI;) = Sat (32) 

together with Eq. ^ implies 

{BiBD^r ^ "^J; (33) 

from which one sees that for each r the Bloch vectors B{,B2, ■ ■ ■ are the vertices 
of a regular d— 1 dimensional simplex. So the d+ 1 families of orthogonal projectors 
define d+1 regular simplices, each having its vertices in BnSo- One might guess, and 
detailed calculation confirms [331 E] i that the optimal choice from a tomographic 
point of view is to choose the projectors in such a way that the polytope formed by 
all d^ +d Bloch vectors has maximal volume. This is achieved if the d+1 different 
simplices are mutually orthogonal: 

{K,Bt)=0 (34) 

for all a, b and r ^ s. 

This condition is often stated in a slightly different form. Suppose we choose 
vectors l-^a) G H such that — li^Dii^al (so for each r the set |V'2)7 • ■ ■ : IV'd) 
is an orthonormal basis for Ti.). Then the requirement that the simplices corre- 
sponding to different bases be mutually orthogonal is equivalent to the requirement 
that 

\{ram\ = ^ (35) 

for all a, b and r ^ s. A family of orthonormal bases for which this condition is 
satisfied is said to be mutually unbiased. 

The question now arises: do families of d + 1 mutually unbiased bases (MUBs) 
actually exist? This is a difficult geometrical problem. As Bengtsson and Erics- 
son [H [Sj have noted, what makes it hard is, in essence, the fact that B D So has 
a much lower dimension than So- Consider, for instance, the case c? = 3. In that 
case the problem is to orientate a set of 4 mutually orthogonal equilateral triangles 
in such a way that all 12 vertices lie in S n So- It is very easy, almost trivial, to 
construct a family of 4 mutually orthogonal equilateral triangles with vertices on 
the 7 dimensional sphere So- The difficult part is then to rotate them so that every 
vertex lies on the 4 dimensional subspace B H So- As it happens the problem has 
been solved for d = 3, and also for every other dimension which is the power of a 
prime number '31', '33J ■ But for values of d which are not prime powers the question 
is stiU open ^di So we have here 

an important physical problem the solution to which depends on gaining a better 
understanding of the geometry of the Bloch body. 

Let us now turn to a different measurement scheme. Suppose that, instead 
of using d+1 different von Neumann measurements, we wanted to use a single 
POVM measurement. The POVM would obviously need to have the property that 
specifying the probability of each of the distinct outcomes fixes the quantum state. 
Such a POVM is said to be informationally complete [Til [T51 [TB I [T7 1 [T ^ [TUl [2U] . 
As we remarked earlier, a complete specification of the quantum state requires the 
specification of — 1 independent numbers. Taking into account the normalisation 
condition (the fact that the probabilities must sum to unity) this means that an 
informationally complete POVM must have at least d^ distinct outcomes. We 
will say that a POVM is minimal informationally complete if it has precisely this 
minimum number of d^ distinct outcomes. The question we have then to answer is: 
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which minimal informationaUy complete POVMs are tomographically optimal? As 
with the MUB problem, the answer to this question depends on achieving a better 
understanding of the geometry of the Bloch body. 

Let El, E2, ■ ■ ■ , E^p be an arbitrary POVM having (P distinct elements. Define 

tr = Tr{Er) (36) 

We may assume that Ej. and consequently for all r (otherwise the POVM 
would effectively reduce to one having fewer than d? elements). It follows that for 
each r the operator (l/tr)Er is a density matrix. We may therefore write, for all r, 

Er^^-^{l+Br) (37) 

where Br G B. The fact that X]r=i ^'r ^ 1 implies 

^tr=d (38) 

r=l 

and 

^trBr=0 (39) 

r=l 

It is easily seen that the POVM is informationaUy complete if and only if the Bloch 
vectors Br span su((i). This in turn will be true if and only if the vectors Br are 
the vertices of a — 1 dimensional simplex (typically an irregular simplex) having 
non-zero volume. One might guess, and detailed calculation confirms [41], that the 
POVM would be optimal from a tomographic point of view if we could arrange 
that (a) the simplex is regular and (b) the vertices all lie on S fl iSo (because the 
volume of the simplex would then be maximal). In other words we would like to 
arrange that 

{Br,Bs) = \'^ , (40) 
' r ^ s 



In that case 



Y,{Br,Bs)^0 (41) 

r=l 

for all s. Since the vectors Bs span sn{d) this means 

'^Br = (42) 

Eqs. ([55]) . ([5^ and (|^^ . taken in conjunction with the fact that the vectors Br 
span the — 1 dimensional space su((i), then imply 

ti=t2 = ■■■ = td2 = ^ (43) 
so that the POVM elements take the form 

Er = -^{l + Br) (44) 

Since the Br all belong to B H So we may alternatively write 

Er = ^Pr (45) 

d 
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where the Pr are a family of one dimensional projectors satisfying 

[d+i ^^'"^ 

The converse is also tru^l: if Pi , -P2 , ■ ■ • , Pd'^ is any family of one dimensional pro- 
jectors satisfying Eq. then ^A, • ■ • j ^-Pd^ is an informationally complete 
POVM. 

A POVM which satisfies the defining Eq. (gHl) (equivalently: a POVM which is 
rank 1 and which satisfies Eqs. and ([^5]) ) is usually referred to as a SIC-POVM 
(symmetric informaationally complete POVM). It appears to us that this terminol- 
ogy is unsatisfactory as, besides being symmetric and informationally complete, 
POVMs of the type in question are also rank 1. As we will see in Section 21 there 
do exist POVMs which are symmetric and informationally complete but not rank 
1. We therefore suggest that POVMs of the type in question would be better de- 
scribed as SI(l)-POVMs ("S" for symmetric, "I" for informationally complete, "1" 
for rank 1). The larger class of POVMs, which are symmetric and informationally 
complete but not necessarily rank 1, we will refer to as SI-POVMs. 

Do SI(l)-POVMs exist? This is a difficult geometrical problem. Moreover, it 
is difhcult for essentially the same reason that the MUB problem is difficult lU [5] : 
namely, the submanifold B (1 So has much lower dimension than the sphere So if 
d > 2. It is easy to consruct a regular cP — l dimensional simplex with vertices in the 
cP ^2 dimensional sphere So, but very hard then to rotate the simplex so that every 
vertex lies in the 2{d—l) dimensional subspace SHiSq (except, of course, when d = 
2). Moreover the difficulty increases with increasing d (because 2{d—l)/{d'^ — 2) 
as d ^ 00). SI(l)-POVMs have been constructed analytically ^Tl [51 [71 [51 ril 
in dimensions 2 to 10 inclusive, and in dimensions 12, 13 and 19. They have been 
constructed numerically [T] in dimensions 5 to 45 inclusive. It is an open question 
whether they exist in dimensions > 45. 

Most (not all) of the SI(l)-POVMs which have been constructed to date are co- 
variant under the action of the Weyl-Heisenberg group (or generalized Pauli group, 
as it is sometimes called). For a summary of the pertinent facts concerning this 
group see Appendix \Xi 

Let be the set of integer pairs p = (^1,^2) such that < pi,P2 < d — 1, 
and for each p S let Dp be the corresponding Weyl-Heisenberg displacement 
operator, as defined by Eq. Let B be any Bloch vector e Bn So- Then the 

fact that the Dp are unitary means that for each p 

Bp = %BDl (47) 

also belongs to BDSo- Suppose that the Bp constitute a regular simplex: 

{Bp,B^) = r ^ P = 'l (48) 

Then the corresponding SI(1)-P0VM is said to be Weyl-Heisenberg covariant. 



■^To see this note that Eq. 1 146 I I implies that the corresponding Bloch vectors satisfy Eq. 14011 . 
It follows that the Bloch vectors span su(d) (because if M is the (d? — 1) X {d? — 1) matrix with 

elements Mrs = (Br,Bs) lor r, s = 1,2, . . . , - 1 then Det M = d^^"^"^ / (cf - l^"^-^ ^ ) and 
,2 

consequently that X]r=i ~ ^ ('-'y same argument that led to Eq. (|42}). The claim is now 
immediate. 
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4. SI-POVMs IN General 

In the last section we discussed SI(l)-POVMs: POVMs which are not only sym- 
metric and informationally complete but also rank-1 (so that each element of the 
POVM is proportional to a one dimensional projector). We now want to broaden 
the discussion, and consider POVMs which, though symmetric and informationally 
complete, are not necessarily rank 1. 

Consider an arbitrary POVM Ei, E2, ■ ■ ■ En defined on a d dimensional Hilbert 
space. Without loss of generality it may be assumed that Er 7^ for all r. We saw 
in the last section that we can write 

Er=^-^il+Br) (49) 

where Br (z B for all r, where i,- > for all r, and where 

J^tr^d (50) 

r 

'^trBr=0 (51) 
r 

Conversely, if we have a set of Bloch vectors Br and positive numbers tr satisfying 
these conditions then Eq. (|49|) defines a POVM. 

We say that the POVM is informationally complete if the probabilities Tr(/5i?r) 
completely specify an arbitrary density matrix p. We say that it is symmetric if 

TiiErEs) = a + f3Sr, (52) 

for all r, s and fixed numbers a, (3. 
We then have the following theorem: 

Theorem 5. Let Ei, E2, ■ ■ ■ , En be a POVM having n elements (all non-zero) 
defined on a d- dimensional Hilbert space. The POVM is symmetric and informa- 
tionally complete if and only if 

(1) n = d^. 

(2) The POVM elements are of the form 

Er = ^{l + Br) (53) 

where the Bloch vectors Br satisfy 

, - ' , I r = s , , 

{Br,B,) = { ^, (54) 
r ^ s 



with < K < 1. 



Remark. We will refer to k as the efficiency parameter as it determines the volume 
of the regular simplex spanned by the Bloch vectors Br, and consequently the 
efficiency of the POVM for tomographic purposes [H]. The POVM is maximally 
efficient if and only if k = 1 in which case it is rank one (an SI(1)-P0VM in the 
terminology explained in the last section). 

Proof. We first prove necessity. Suppose the POVM is symmetric and information- 
ally complete. We can write it in the form 

Er^^-^{\+Br) (55) 

a 
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for Bloch vectors Br and positive numbers tj. satisfying Eqs. ([50]) and ((5T|) . The 
symmetry condition Eq. (|52p then imphes 



TT:{Er) ^ T^r{ErEs) = na + f3 (56) 



s=l 



for aU r. In view of Eq. (|50p this means 



= - (57) 
n 



(58) 



for all r, and consequently 
Using these results, Eq. and the symmetry condition Eq. (|52p we deduce 

The fact that the Br are Bloch vectors means {Br, Br) < 1. We must also have 
{Br, Br) > (because otherwise Er = for all r, in which case the POVM would 
not be informationally complete). Consequently 

0</3<^^ (60) 
n(n — 1) 



Let M be the nx n matrix with elements Mrs = {Br,Bs). Since the POVM is 
informationally complete the Bloch vectors Br must span the cP — 1 dimensional 
space su(c?). So M must have rank d? — 1. On the other hand 

Det(M-A)^-A(-^-A) (61) 

It follows from this that M has n — I non-zero eigenvalues (since we have shown 
that /3 > 0). However, the fact that M is rank cf — 1 means that it must have 
(P ~ 1 non-zero eigenvalues. We conclude that n = cP. Making the substitutions 
n — (P and /3 = ^^^^^^-^ in Eq. ([S^ we obtain Eq. (|54p . Moreover, it follows from 
Eq. (Uni) that < k < 1. 

Having proved necessity, it remains to prove sufficiency. Suppose Bi, B2, ■ ■ ■ , B(i2 
are Bloch vectors satisfying Eq. (IM)) . Let M be the <P x (P matrix with elements 
Mrs = {Br,Bs). Then 

Det(Af-A) = -A(-^^-Aj (62) 

Since, by assumption, k > it follows that M has (P — 1 non-zero eigenvalues, and 
is therefore rank (P~l. Consequently the Bloch vectors span the (P — 1 dimensional 
vector space su((i). 
Eq. also implies 

'(j2Bs),Br)^0 (63) 

s=l 

for all r. Since the Br span sn{d) we deduce 



Y^Bs^O (64) 
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It follows from this that if we define 

1 

d2 



Er^—{1 + Br) (65) 



the operators Ei, E2, . . . , Ed^ constitute a POVM. The fact that the Br span su{d) 
means the POVM is informationally complete. The fact that the POVM is sym- 
metric is immediate. □ 

We noted in the last section that the existence problem for SI(l)-POVMs is 
hard, and still unsolved for dimensions > 45. But if one relaxes the demand that 
the POVM be rank 1, and simply looks for an SI-POVM of arbitrary rank, the 
problem becomes much easier. 

To construct an SI-POVM of arbitrary rank all we have to do is construct a 
regular simplex in su(d) with its vertices all on So (since So is a sphere such simplices 
are guaranteed to exist). Let Bi,B2, - ■ ■ B^p be the vertices. Then 

{Br,Bs) = \^ , '\[^'. (66) 

I — otherwise 

If the Br were Bloch vectors this would give us an SI(1)-P0VM. However, if the 
simplex is chosen at random they are very unlikely to be Bloch vectors (because the 
manifold Bf^So has much lower dimension than So)- Nevertheless, we can still use 
them to construct an SI-POVM by shrinking the simplex until the vertices are all 
in B. In fact, let —rrir- be the smallest eigenvalue of Br- It follows from Thcorcm[l] 
that 1 < rrir- < d — 1 for all r. Now define 



K = mm 

l<r<d2 



(^) 



We have -^-^ < « < 1- Moreover, it follows from Theorem [2] that B'r = nBr G B 
for all r. By construction 

(B'r,B':) = \ . (68) 

I— ^puj- otherwise 

So we can use Theorem [5] to deduce that the POVM with elements 

Er^^{\^B'r) (69) 

is symmetric, informationally complete with efhciency parameter = k. 

The argument just given shows that in every dimension d there exists an SI- 
POVM with efhciency parameter > ^j^. We will see in Section [7] that at least 
when d is odd it is possible to considerably improve on that. 



5. SI-POVMS WHICH ARE Weyl-Heisenberg Covariant 

In Section [6] we will discuss the bearing of the above results on the really difficult 
problem: i.e. the problem of constructing POVMs which are, not merely symmetric 
and informationally complete, but also rank 1 (have efficiency parameter = 1). In 
preparation for that we first need to prove a result concerning SI-POVMs (with 
efficiency parameter not necessarily — 1) which are covariant under the Weyl- 
Heisenberg group. 

We begin with a definition. Let B E So (we do not assume that B is a Bloch 
vector), and for each p G let Bp = DpBD^ (where and -Dp are as defined 
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in Appendix |X| . We say that B is the generating vector for a Weyl-Heisenberg 
covariant regular simplex if 

, ''-^<°-»' (70) 

I — ^|33j otherwise 

It is easily seen that if that is the case 

I — otherwise 

meaning that the vectors Bp are the vertices of a regular simplex. 
We now have the following lemma: 

Lemma 6. A vector B d So is the generating vector for a Weyl-Heisenberg covari- 
ant regular simplex if and only if 

^^Vm ^ (72) 

for any set of real numbers 6q satisfying the condition e**^ = S-^e^^^'^ (where (Z^)*, 
s_q and q are as defined in Appendix C^) . 

Proof. We know from Eq. (|112p that any vector B E So can be written 

where the expansion coefficients Cq = {l/d)Tr{D^B) satisfy the condition Cq = 
s_qC* . By a straightforward application of Eq. (j95p we find 

In view of Lemma [7] in the Appendix it follows 



Suppose now that |cq| = l/\/d + 1 for all non-zero q. Then Eq. fTS]) implies 

(5,Sp> - ^(-1+ E = ^,i~^ + dXo) (76) 

So i? is the generating vector for a Weyl-Heisenberg covariant regular simplex. 
To prove necessity, suppose that Eq. ((75)) is satisfied. Using the fact that 



for all q, r G to invert Eq. ([TS]) one finds |cq| = l/Vd+T for all non-zero q. □ 

This lemma gives us an easy way to construct SI-POVMs. Simply choose an 
arbitrary set of phases e*^'^ satisfying the condition e*^"! = s_qe~**'' and construct 
the vector B specified by Eq. ([72]) . Let — 1/k be the minimum eigenvalue of B. It 
follows from Theorem[T]that 1/ (d—l) < k < 1. Moreover —1/k is also the minimum 
eigenvalue of Bp for all p. So it follows from Theorem [5] that the operators 

4 = ^(l + «Bp) (78) 
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consitute a POVM. By construction the POVM is SI, Weyl-Heisenberg covariant, 
and has efficiency parameter k > \/{d— 1). 



6. Construction of SI(1)-POVMs 

Of course, what we would really like to do is to construct a POVM which is, 
not merely symmetric and informationally complete, but also rank 1. The POVM 
defined by Eq. (|78p will be rank 1, with efficiency parameter k = 1, if and only if 
5 is a Bloch vector. The question therefore arises: how do we choose the phases in 
Eq. fT^]) so as to ensure that that is the case? 

We can answer that question by appealing to Theorem [H The vector B in 
Eq. ([7^ is on the sphere Sq. So Theorem H] tells us that 

Tt:{B^) < d{d - l){d - 2) (79) 

with equality if and only if _B is a Bloch vector. In terms of the phases on the right 
hand side of Eq. ([7^ the condition reads (using Lemma [7] in the Appendix) 

E Sp+qr<P''i>e'(''p+''''-''-®<') < (d - l){d ~2){d+ 1)^ (80) 

p,q,peqG(Z2)' 

with equality if and only if i? is a Bloch vector. 

This gives us an extremal condition alternative to the one used by Renes et al [T] . 
Renes et al [T] base their numerical construction of Weyl-Heisenberg covariant SI(1)- 
POVMs on the fact that, if P is an arbitrary rank 1 projector and Pp = DpPD^, 
then 

E (Tr(^pP))' > ^ (81) 
pez^ 

with equality if and only if the operators jPp constitute an SI(1)-P0VM. The 
inequality we have derived provides us with an alternative procedure: instead of 
looking for a projector P which minimizes the expression on the left hand side of 
Eq. (|8ip. one can look for a set of phases which maximize the expression on the left 
hand side of Eq. 

It should be said that if one is specifically looking for a method of constructing 
SI(l)-POVMs numerically a procedure based on Eq. ([5T|) is likely to be more efficient 
than one based on Eq. ([50]) . This is because the expression on the left hand side 
of Eq. ([81]) is a function oi2{d— 1) real parameters {i.e. the number of parameters 
needed to specify the projector P), whereas the one on the left hand side of Eq. (150]) 
is a function of (d^ + 1)/2 real parameters if d is odd and (d^ + 4) /2 real parameters 
if d is even (i.e. the number of independent phase angles). 

However, although the extremal condition represented by Eq. ([501 would appear 
not to have any advantages from a concrete numerical point of view, it may per- 
haps be interesting from a more abstract mathematical point of view, as providing 
additional insight into the problem. In particular, the fact that the phase angles 
appear in combinations of the form Op + dq + Or with p -I- q -I- r = (mod d) may 
possibly provide some clues as to the origin of the order 3 symmetry found in every 
Weyl-Heisenberg covariant SI(1)-P0VM constructed to datc0. 



^Note that Grassl [8] has constructed a counter-example in dimension 12 to conjecture C of 
ref. [7]. However, his example is still invariant under a canonical order 3 unitary. Specifically 



his matrix T12 is a representative of the Clifford operation 



which it will be 



seen has Clifford trace = — 1 (notation and terminology as in ref. ^7^). His example is therefore 
consistent with conjecture A of ref. [?]• 
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7. The Wigner POVM 

Suppose that d is odd. In that case we can set the phase angles on the right 
hand side of Eq. ([72| equal to zero, giving 

^ = ^TTf E ^1 (82) 

(notice that if d was even this choice of phases would not be permissible because 
when d is even the signs s_q are not all positive). For reasons explained below we 
will refer to the SI-POVM corresponding to this choice of B as the Wigner POVM. 

We wish to determine the efficiency parameter of the Wigner POVM. For that 
purpose it is convenient to consider the operatoiQ 

U =hl + VdTlB) = ^y (83) 
d d 

U, like is an Hermitian operator. Moreover 

= 1 (84) 
where we used the fact that X^rez^ t^'^'''^ = d^6qo (note that this depends on the 

fact that d is odd). It follows that the eigenvalues of [/ all = ±1. Taking into 
account the fact that Tr(C/) = 1 we deduce that U must have {d+ l)/2 eigenvalues 
= 1 and {d — l)/2 eigenvalues = — 1. Consequently the smallest eigenvalue of B 
is —\/d + 1. In view of Theorem [2] it follows that {1/^/d + 1)B is a Bloch vector. 
Hence the d^ operators 

constitute an SI-POVM of rank (d+l)/2. We will refer to this as the Wigner POVM. 
It has efficiency parameter 1/vd+l — which is a considerable improvement on the 
worst case value l/{d— 1) calculated in Section [51 although still greatly inferior to 
the best case value k = 1. 

Let us now explain the connection between the Wigner POVM and the Wigner 
function. Let p be an arbitrary density matrix, and let 

Pp = i Tr{Dlp) (86) 

We define the Wigner function Wp to be the discrete Fourier transform of the 
coefficients pp: 

Wp^^Y. ^''^"^""^Pp (87) 

This definition agrees with that of Wootters [52] in the case when d is prime. If 
d is non-prime the Wigner function as defined by this formula loses some of the 
properties which Wootters considers desirable. However, it appears to us that it 



'^In the notation of ref. ^ U is a representative of the CUfford operatic 



-1 0\ 0-) 

Its action on the standard basis used to define the operators Dp (see Eqs. II91I I and II92I I) is 
U\r) = |f}. So IJ can be thought of as a discrete parity operator. 



15 



retains sufficiently many of these properties for it still to be considered a reasonable 
way of defining the Wigner function. For further discussion of the discrete Wigner 
function see refs. [32l [Ml HSl HSl HH |45] and references cited therein. 

The Wigner function can be expressed in terms of the operators Up = DpUD^ 
(where U is the operator defined in Eq. ([S^]) ). In fact 

Eqs. (IMl) and ([87]) then imply 

Wp = ^ Ti{Upp) (89) 

Taking into account Eqs. ([83]) and ([85]) we deduce 

Wp = {d+l)TT(Epp)^^ (90) 

Of course, the fact that the Wigner function is a linear function of the probabilities 
Ti{Epp) is an automatic consequence of the fact that the POVM is informationally 
complete. However, in the case of the Wigner POVM the relationship is particularly 
simple: to obtain the Wigner function one merely has to rescale the probabilities 
by a constant amount and then shift them by a constant amount. 

Eq. ([^ is conceptually interesting because it establishes a connection between 
SI-POVMs and the Wigner function. At first sight it may appear that it also 
has a more concrete, pragmatic significance, as providing a good way to deter- 
mine the Wigner function tomographically. However, a little reflection will dispel 
that impression. The trouble is that the Wigner POVM has efficiency parameter 
= l/VcT+T, which is < 1 (and ^ 1 if d is large). So if one wants to determine the 
numbers Wp it would be much more efficient (would give much less statistical un- 
certainty for a given number of measurements) to use a tomographic scheme based 
on an SI(1)-P0VM, or a full set of MUBs (in dimensions where such exist), and 
then to perform the appropriate linear transformation on the relative frequencies 
obtained [H]. 

Finally, let us note that Miquel et al ^42j have described a scheme for "directly 
measuring" the individual numbers Wp. This scheme might, perhaps, have some 
advantages over a scheme based on an SI(1)-P0VM or a full set of MUBs in a case 
where one was only interested in some of the numbers Wp. 

8. Conclusion 

We originally undertook the investigations reported here in the hope that they 
might lead to a solution of the really challenging problem, which is to demonstrate 
the existence (or, as it may be, the non-existence) of SI(l)-POVMs in every finite 
dimension. We did not succeed in that primary aim. Nevertheless, we derive some 
consolation from the fact that the class of SI-POVMs is of some intrinsic interest. 
Also, it is not impossible that the results reported here contain clues that may help 
us to solve the main problem. 

Appendix A. Weyl-Heisenberg Group 

In this appendix we summarise those facts concerning the Weyl-Heisenberg group 
(or generalized Pauli group as it is sometimes called) which are needed in the main 
text. Our definitions are those of ref. [3, and may differ slightly from the ones used 
by other authors. Let |0), . . . jd — 1) be an orthonormal basis for Ti., and le10 

^The reason for defining r = —e^^l'^ rather than r = e'^*/'' is that it means t'^ =1 for all d. 
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T = —e^^l'^. Define operators T and S by 

T\r) = r^^r) (91) 

MO) r = a — 1 

Then define, for each pair of integers p = (pi,P2) G 

Dp ^ tP'P^ SP^fP^ (93) 

Tfie operators -Dp are the displacement operators of the Weyl-Heisenberg group. 
The reason for including the factor tP^p^ is that it means that the operators have 
the following nice properties: 

Dl = D_p (94) 
DpD^^T<P'^^Dp+^ (95) 

where 

(p, q) = P2qi - piq2 (96) 

The fact that (p, p) = means 

(^p)" = -D„p (97) 
for all p € and n € Z. In particular the operators Dp are unitary: 

DlDp = 1 (98) 
for all p. It is also worth noting that 

(^p)' = l (99) 
for all p (this is one of the reasons for setting r = — e*^/'*. If, instead, one set 
T = e"/'' it would sometimes happen that (-Dp)'' = ^1)- 

The presence of the factor t^P''!^ on the right hand side of Eq. (PSj) means that 
the operators -Dp do not constitute a group. However, one obtains a group (the 
Weyl Heisenberg group) if one takes the set of all operators of the form e*"-Dp, 
where e*" is an arbitrary phase (alternatively, one can define the Weyl-Heisenberg 
group to be the set of all operators of the form r^Dp, where n is an arbitrary 
integer). 

If p = q (mod d) then .Dp = Dq up to a sign. Specifically: 

I)p = P\., ^^'^^^^'^^ (100) 

|(-l)3<P^i>-Dq if d is even ^ ' 

(to prove this formula write p = q + du and then use Eq. (|95p ). It is therefore often 
convenient to restrict ourselves to values of p lying in the set = { (pi j P2 ) : Pi,P2 = 
0, 1, 2, . . . , c? — 1}. Given arbitrary p e Z^ let [p] be the unique element of Z^ such 
that [p] = p mod d. It is also convenient to define 

p®q=[p + q] (101) 

p e q = [p - q] (102) 
P = [-P] (103) 



and 




if d is odd ^^^^^ 
if d is even 
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We then have, for aU p, q G Z^, 

Dl = s^pDp (105) 

DpD^ = Sp+qT<P^1>7)peq (106) 

It is also easily verified that 

Ti{DlD^) - d5pq (107) 

for all p, q G Z^. This means that, relative to the Hilbert-Schmidt inner product, 
the operators "^^p are an orthonormal basis for the cP complex dimensional space 

jC(?^) consisting of all d x rf complex matrices. So an arbitrary matrix A e >C(H) 
can be expanded 

where the expansion coefhcients are given by 

Ap = iTr(l)ti) (109) 

It follows from Eqs. (|105p . pOSp and (|109p that A is Hermitian if and only if 

Ap = s_pA; (110) 

for all p e Z2. 

Let {ZD* = {p e Z2 : p ^ (0, 0)}. The fact that 

Tr(i)p) = |'^ ifp = 0(modd) ^^^^^ 
I otherwise 

means that A e su{d) if and only if it has an expansion 

A^ J2 ^P^P (112) 

pe(z5)- 

where the coefhcients satisfy Eq. (jllOp . 

The following lemma tells us how to calculate the expansion coefficients and 
traces of double and triple products: 

Lemma 7. Let A,B,C e C{H). Then 

(iS)p = ^ Sp^qT^'l-P^^qiSpeq (US) 

{ABC)^= J2 ■'^P-q-rr<'i+"^^P>+<'i'"->Aqi3rCpeqer (114) 

where {AB)p = (1/d) Tr(i)|,ii?) anrf {ABC)p = (1/d) Tr(I)j,ii3C') are f/ie expan- 
sion coefficients as given by Eq. m09\) . Traces are given by 



Tt{AB) = d ^ S-q^qSq (115) 

qGZ^ 

TT{ABC)=d s_q_rr<i^'-UqBrCqef (116) 



q,re2 
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If A,B,C are Hermitian we can alternatively write 

Tr(iS) = d A^B*^ (117) 

TT{ABC)^d J2 Sq+rr<'i^'->AqBrC;er (US) 
Proof. Let A, B e C{n). It follows from Eq. (fTUSl) that 

q,reZ^ 

= ^ Sq+rT^^l'-^Uq^r^qer 

q,r6Z5 

= E Sq+peqT<'»'P®'>>AqBpeq^p 

= E ( E Sp^qr<1 P^qSpeq ) ^p (119) 

where in the last line we used the fact that Sq+peq'''^'''^®'*^ = Sp-qT^'^'P^ . Eq. (|113p 
is now immediate. 

To prove Eq. (I114p we apply Eq. pi3p twice: 

(iBC')p = J2 ^P-^r^^'^^MBQpe^ 

qezi 

= E Sp-qSpeq-rT^'i'P^+^'-'Pe'i^AqSrCpeqer 

= J2 Sp-q-rr<'i+'-'P>+<'i''-Uqi3rCpeqer (120) 



q.reZ^ 

where in the last line we used the identity Sp_qSpQq_rT^'"'P®i^ = Sp_q_r'''^'^'^~''^ • 
To prove Eqs. ITTSl) and (fTTB]) set p = in Eqs. ([11^ and (fTTIl) and use the fact 
that Tt{M) = dMo for all M e C{H). Eq. (fTTTll follows from Eq. (fTTS]) . Eq. (fTTS]) 
follows from Eq. pi6p and the identity Sp+qS-p-q = s_p0q. □ 
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